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SECTION  1 
INTRODUCTION 


The  theory  of  the  electromagnetic  pulse  (EMP)  from  high-altitude 
nuclear  bursts  was  first  developed  in  a series  of  lectures  (Reference  1) 
given  by  this  author  at  the  Air  Fcrce  Weapons  Laboratory  (AFWL)  in  January 
and  February  1964.  The  essential  part  of  the  theory  is  that  Compton  recoil 
electrons,  produced  by  the  prompt  gamma  rays  from  the  burst,  are  deflected 
by  the  geomagnetic  field  from  the  radial  direction  to  a direction  perpen- 
dicular to  both  the  radial  vector  and  the  geomagnetic  field.  The  transverse 
electric  current  so  formed  generates  an  outgoing  EM  wave  which  maintains 
coincidence  with  the  gamma  pulse  and  (therefore)  with  the  Compton  current; 
as  a result  of  this  coincidence,  a large- amplitude  pulse  is  built  up,  with 
duration  (several  tens  of  nanoseconds)  determined  by  Compton  electron 
dynamics.  This  short  duration  is  in  marked  contrast  to  that  of  the  EMP 
radiated  by  the  radial  Compton  current,  which  is  determined  by  the  size  of 
the  source  region  (hundreds  of  kilometers)  and  leads  to  a smaller  amplitude. 


In  Reference  1 the  author  derived  the  basic  equation  of  the  out- 
going wave  approximation. 
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Here  r is  the  radial  coordinate  (but  with  the  gamma  pulse  treated  as  planar), 
is  the  transverse  electric  field,  Jt  is  the  transverse  Compton  current 
density,  a is  the  electrical  conductivity  induced  in  the  air,  and  ZQ  * 377 
ohms  is  the  impedance  of  space.  The  factors  i/2  occur  because  only  outgoing 


(and  not  ingoing)  waves  are  generated  with  substantial  amplitude.  The 
equation  says  that,  as  r increases,  L^  builds  up  in  the  direction  of  -J  , 
but  is  attenuated  by  the  conductivity.  The  retarded  time  t-r/c  enters 
liquation  0 only  as  a parameter;  that  is,  the  equation  deals  with  F.  and  J 
at  a constant  retarded  time.  Sample  Compton  currents  and  conductivities 
were  calculated  in  Reference  1,  and  the  solution  presented.  A general 
discussion  of  the  character  of  solutions  was  given,  along  with  a discussion 
of  the  effect  of  diffraction,  which  explained  why  the  solution  of  an 
equation  along  a single  ray  gives  correct  answers  for  a three-dimensional 
problem. 


In  largely  independent  work  only  slightly  later  than  this  author’s, 
William  Karzas  and  Richard  Latter  developed  virtually  the  same  theory  (Reference 
2).  These  authors  chose  to  work  in  spherical  coordinates,  in  which  the 
derivative  3Et/3r  in  F.quation  0 is  replaced  by  3(rF-t)/r3r,  and  they  developed 
approximate  analytical  formulae  for  the  Compton  current,  the  conductivity, 
and  the  resulting  field.  There  has  been  no  important  disagreement  over  the 
basic  theory  of  the  high-altitude  BMP  since  1964. 

For  a couple  of  years  after  the  birth  of  the  theory,  Karzas  and  this 
author  provided  analytically-based  estimates  of  the  high-altitude  EMP  to 
military  systems  planners.  Then  John  Erkkila  of  AFWL,  working  with  consulta- 
tion from  Karzas,  constructed  the  computer  code  HEMP  which  essentially 
mechanized  the  analytical  models  developed  by  Karzas  and  Latter.  This  code 
made  it  possible  to  produce  EMP  environment  information  easily  for  many  burst 
heights,  yields,  and  observer  locations,  and  it  was  the  principle  source  of 
such  information  for  several  years. 

The  analytical  models  contained  quite  a few  approximations.  In 
1971,  Jerry  Longley  and  this  author  undertook,  with  support  of  the  Defense 
Nuclear  Agency  (DNA),  to  develop  a computer  code  with  more  accurate  modeling 
of  the  basic  phenomena . First,  by  solving  the  equations  of  motion  of  sample 
Compton  electrons  simultaneously  with  the  solution  of  Maxwell ' s equations. 
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we  obtained  (Reference  3)  the  effect  of  the  KMP  fields  on  the  Compton  current 
(self  consistency).  Second,  wc  developed  a method,  called  the  obliquity 
factor  method,  for  treating  the  effect  of  multiple  scattering  of  Compton 
electrons  by  air  atoms  (Reference  4).  Third,  we  developed  a method  for  treat- 
ing the  finite  time  required  for  secondary  electrons  to  complete  their 
i on i cat  ion  (Reference  5).  The  code  which  embodies  these  improvements,  and 
others,  is  called  CHAP  (Reference  6).  The  AFKc  codes  CllliMP  and  lll-MP-B 
employ  the  same  methods,  and  most  of  the  current  high-altitude  BIP  environ- 
ments are  computed  with  one  or  another  of  these  codes. 

Recently  (Reference  7),  William  Sollfrcy  has  calculated  the  effect 
of  multiple  scattering  on  the  Compton  current  by  a new  method  developed  by 
him,  and  has  raised  questions  about  the  basis  and  the  accuracy  of  the 
obliquity  factor  method. 

In  the  present  report  we  explain  the  basis  of  the  obliquity  factor 
method,  and  test  its  accuracy  against  Monte  Carlo  calculations.  We  shall 
see  that  it  is  quite  accurate.  It  is,  in  fact,  more  accurate  than  we  had 
expected. 


We  also  «'ompare  with  Sollfrey’s  calculations  and  find,  to  our  dis- 
appointment, that  his  method  does  not  apparently  provide  accurate  answers 
for  this  difficult  problem  of  multiple  scattering.  We  are  therefore  left 
with  having  to  rely  on  the  Monte  Carlo  calculations  for  accuracy  standards. 


SECTION  2 

PROBLEM  DEFINITION 


The  prompt  gamma  rays  from  a nuclear  burst  are  emitted  in  a few 
tens  of  nanoseconds.  For  a burst  far  above  the  atmosphere,  the  gamma  rays 
at  time  t after  the  burst  lie  within  a spherical  shell  of  radius  ct  (c  is  the 
speed  of  light)  and  with  thickness  of  the  order  of  10  meters.  The  downward 
going  part  of  this  shell  begins  to  interact  with  the  atmosphere  at  an 
altitude  of  about  50  km.  By  30  km,  where  the  gamma  scattering  length  is  about 
equal  to  the  atmospheric  scale  height  h 6.7  km,  of  the  order  of  one-half 

of  the  gammas  have  been  scattered  by  the  Compton  scattering  process.  By 
20  km,  only  about  1 percent  of  the  gammas  have  not  been  scattered.  The 
dominant  source  region  for  the  high-altitude  liMP  lies  between  20  and  40  km 
altitude. 


The  mean  scattering  angle  0 of  the  gammas  is  of  the  order  of 
30  degrees.  Hence  in  the  interval  before  a second  scattering  of  a gamma 
occurs  it  will  fall  behind  the  uascattered  gammas  by  a distance  ~ A(l-cosO)  t* 
0.13  1 k l km.  The  scattered  ga.uui-.  therefore  occupy  a much  thicker  shell 
than  the  unscattered  gammas,  and  tt.o  flux  of  scattered  gammas  is  very  small 
compared  with  that  of  the  unscattered  gammas  in  the  thin  shell  occupied  by 
the  latter.  For  the  first  hundred  nanoseconds  of  the  F-MP,  the  previously 
scattered  gammas  can  be  ignored. 

The  mot  ion  of  Compton  recoil  electrons  is  limited  by  their  Farmer 
radius  in  the  geomagnetic  field,  and  by  energy  loss  in  the  air,  to  distances 
of  the  order  of  100  meters  fro*  their  birth  place.  Over  this  distance  one 
can,  with  little  error,  consider  the  gamma  shell  to  be  planar,  the  gammas 
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parallel,  and  th  r intensity  constant  in  space  at  a given  retarded  time 
I , 

I = t - 2 . ( I ) 

c 

Ik-rc  he  have  ust-d  the  Cartesian  coordinate  z instead  of  r as  t lie  < large) 
distance  from  the  Imrst  point.  (Of  course,  once  we  have  gone  to  Cartesian 
coordinates,  we  can  choose  tin-  origin  of  z arbitrarily.)  Thus  over  tin- 
region  of  space  of  interest  in  calculating  the  foiupton  current,  we  can  regard 
the  gamma  flux  l\  as  being  a function  of  T alone, 

I = F IT)  . (2) 

1 ( 

Since  the  air  density  and  geomagnetic  field  are  also  very  nearly  constant 
over  100  meters,  the  Compton  current  density,  the  air  conductivity,  and  the 
IMP  fields  are  also,  to  good  approximation,  functions  of  the  retarded  time 
alone  over  the  range  of  the  Compton  electrons.  This  approx imat ion  greatly 
simplifies  the  problem  of  calculating  the  Compton  current  and  the  F.MP. 

The  EHP  fields  affect  the  motion  of  the  Compton  elect rons.  Thus 
for  an  accurate  calculation  it  is  necessary  to  solve  Maxwell's  equations  for 
the  fields  simultaneously  with  the  equations  of  motion  of  the  Compton  electrons. 
Since  the  combined  problem  of  particles  and  fields  is  nonlinear,  there  is  little 
hope  of  solving  it  accurately  by  analytical  means  for  general  cases;  we  are 
reduced  to  using  numerical  methods. 

In  considering  numerical  methods,  we  have  to  choose  between  dealing 
with  particles  or  with  a distribution  fimction,  i.e.,  between  solving  Newton's 
laws  or  the  Boltzmann  equation.  In  the  Boltzmann  equation  treatment,  the 
distribution  function  will  be  a function  of  four  variables:  T and  the  three 
components  of  electron  momentum.  The  fine  gridding  of  momentum  space  re- 
quired for  accurate  solutions  makes  this  method  generally  more  time  consuming 
thar.  particle  methods.  (Ne  do  not  want  to  rule  out  the  possibility  of 
inventing  a clever  method  for  reducing  the  number  of  grid  points  required.) 
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we  tii ink  lirst  of  Monte  farlo 


i. lining  to  particle  t reatment s , 
met  hods , because  of  tli<-  phenomenon  of  scat  tor i n j»  of  the-  fompton  electrons 
|,y  air  atoms,  which  must  Iji:  included  for  an  accurate  solution.  However, 
the  random  fluctuations  in  Monte  f.arlo  calculations  will  falsify  the  f* i j»l» - 
frequency  content  of  the  calculated  IMP  unless  very  large  numhers  of  electrons 
are  piucessed . "I  he  computing  time  consumed  males  the  Monte  f.arlo  method 
suitable  only  for  infrequent  cheers  on  faster  methods. 

It  can  lie  estimated,  and  we  shall  later  show,  that  scattering 
causes  only  about  a Mi  percent  reduction  of  the  peak  Compton  current  at 
Mi  km  altitude,  which  is  the  center  of  the  IMP  source  region  and  the  place 
where  the  largest  IMP  is  generated . Thus  there  is  hope  that  an  approximate 
method  of  including  scattering  could  give  answers  accurate  to,  say,  10  percent, 
winch  would  he  adequate.  To  this  end  we  invented  the  obliquity  factor 
method.  It  includes  the  effect  of  election  scattering  without  use  of  random 
numbers . 
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SECTION  3 

COMPTON  SCATTERING 


The  Compton  recoil  electrons  (Reference  8)  are  produced  in 
collisions  of  the  gamma  rays  with  electrons  in  the  air  atoms.  In  such  a 
collision,  the  gamma  is  scattered  to  angle  0 from  its  original  direction 
(see  figure  1),  and  the  electron  goes  off  at  angle  i/i.  The  initial  and 
final  gamma  directions  and  t he  electron  direction  lie  in  a plane.  The  azimuthal 
angle  of  this  plane,  or,  to  he  specific,  of  the  scattered  gammas,  about  the 
original  gamma  direction  is  denoted  by  <j>.  We  neglect  the  effect  of  binding  of 
the  electrons  in  the  atoms,  which  is  believed  to  be  small. 

We  shall  use  the  convenient  relativistic  notation  in  which 
•> 

Y = photon  cnergy/mc"  = photon  momentum/me 

f;  = electron  total  energy/mc*-  l (3) 

p = electron  momentum/mc  ' 


Figure  1 . Geometry  of  Compton  scattering. 
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Mere  m is  the  electron  rest  mass,  and  c the  speed  of  light.  From  energy 
and  momentum  conservation  one  can  determine  the  dependence  of  y ' and  f,  oi 
y and  0: 


1 + 7X  ’ 


i = 1 + 7 - y'  = 1 + ■ , 

1 + yy  ’ 


where 


X = 1 ~ cosB  , 0 £ y<  2 . 

One  also  finds  the  z-component  of  electron  momentum 


(4) 

(5) 


(6) 


pz  = y - y'cose  = y(i+y)  y-jh^  • (?) 

The  additional  relativistic  relations, 

c = V P2  ♦ 1 , (S) 


v 


_z 

c 


e 


etc.  , 


(9) 


are  also  useful.  In  Equation  9,  vz  is  the  z-component  of  electron  velocity; 
similar  equations  hold  for  the  other  components. 


The  probability  of  scattering  to  angle  0 is  given  by  the  formula 
of  Klein  and  Nishina.  The  differential  cross  section  of  an  electron  for 
scattering  the  gamma  into  solid  angle  dft  at  0 is 


o(Q)dQ  = 


dfi  [\  „ ,2 

5-  1 + (1-x)  ♦ 

Unx)  l 


__A._ 

1 ♦ YX 


(10) 


2 2 

Here  rQ  = e /me  is  the  classical  electron  radius,  where  -e  is  the  electron 
charge  (cgs  units).  The  element  of  solid  angle  is,  after  integrating  over 
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d(i  = 2trsin0d0  = 27id/  . 


HI) 


The  total  scattering  cross  section  per  electron  is  obtained  by  integrating 
liquation  10  over  x>  and  is 

nr  = tirjyi'jfYj  , (12) 

where 


T,(Y) 


2(2+«y+9y2+T3)  2 + 2y  - / 

2 '' 2 ’ " t 

Y f 1*2Y3  Y 


9m  ( l + 2y) 


lor  the  gamma  flux  given  by  liquation  2,  the  total  source  density 
of  Compton  electrons  is 

St(T)  = NZotI-y(T)  . (14) 

Here  N is  the  density  of  air  atoms  and  Z = 7.2  is  the  mean  atomic  number  of 
air. 
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SECTION  4 

THE  BOLTZMANN  EQUATION 


We  originally  derived  the  obliquity  factor  method  by  considering 
particles  ru  er  than  the  distribution  function.  Since  that  derivation 
apparently  left  something  to  be  desired  in  clarity,  we  shall  derive  it  here 
from  the  Boltzmann  equation.  Let  f(r,p,t)  be  the  density  of  L'ompton 
electrons  in  the  six  dimensional  phase  space  of  coordinate  r and  momentum  p. 

The  six-dimensional  velocity  of  the  particles  in  phase  space  near  the  point 
r,p  is  r,p,  where  the  dots  indicate  time  derivatives  of  the  particle 
quantities: 

r = v = three  dimensional  particle  velocity,  (15) 

p = force  on  particle  at  r,v  . (16) 

Part  of  the  force  comes  from  the  electric  and  magnetic  fields.  We  shall 
also  imagine  that  a resistive  force  acts  on  the  Compton  electrons,  to  account 
for  their  gradual  loss  of  energy  to  other  electrons  in  air  atoms.  The 
scattering  of  the  Compton  electrons  by  air  atoms  will  be  treated  as  a 
separate,  scattering  term  in  the  Boltzmann  equation. 

The  six-dimensional  current  of  particles  in  phase  space  is  (r,p)f, 
and  the  six-dimensional  gradient  operator  is  (7  ,V  ) where  the  subscripts 

* r 

indicate  whether  the  three-dimensional  gradient  operates  on  coordinate  or 
momentum  space  variables.  The  six-dimensional  divergence  of  the  six-dimensional 
particle  current  is 


(17) 


CVV  ' [(?'P)f]  = vr  ■ (7f)  + Vp  ■ fpf) 

= v • Vrf  * vp  . (pf)  . 

Here  the  order  of  and  v can  be  interchanged  because  v is  a function  of 
the  momentum  variables  (see  liquations  8 and  9),  which  are  independent 
variables  from  r.  The  order  of  V and  p cannot  be  interchanged  because  of 
the  resistive  force  which  depends  on  p in  such  a way  that  V • p ^ 0. 

The  conservation  of  particles  is  expressed  by  the  Boltzmann 

equat ion, 

||  + v • Vrf  ♦ 7p  • (pf)  = S + fK«’„CP’)dV  • (18) 

Here  S is  the  source  density  of  Compton  electrons  in  Compton  collisions, 
and  the  integral  with  K is  the  scattering  operator,  which  takes  particles 
out  of  momentum  p'  and  places  them  at  momentum  p. 

It  is  important  to  realize  that  the  Boltzmann  equation  is  completely 
equivalent  to  Newton's  laws  of  motion.  If  we  start  with  an  f which  is  a sum 
of  delta  functions,  each  singularity  representing  a point  particle,  and  solve 
the  equation  forward  in  time,  then  the  delta  functions  will  be  preserved 
and  they  will  move  exactly  as  particles  would  under  Newton's  laws.  To  make 
this  work,  the  source  S and  the  scattering  operator  have  to  be  regarded  as 
stochastic  operators,  which  occassional ly  inject  additional  point  particles 
or  scatter  point  particles  from  one  momentum  to  another.  On  the  other 
hand,  one  can  regard  f as  a continuous  function  expressing  the  probability 
of  finding  a particle  near  r,p,  or  as  the  density  of  particles. 

We  have  seen  that  it  is  a good  approximation  to  treat  the  source 
S as  depending  on  r and  t only  through  the  retarded  time. 


(19) 


S = S(t  -|,p)  = S(T,p)  . 

With  this  source,  the  Boltzmann  equation  and  Maxwell's  equations  allow 
solutions  which  are  functions  of  T alone.  On  substituting 

f = f(t  ,p)  , (20) 

into  the  Boltzmann  equation  we  find 

(1  -~)f  ♦ • (pf)  = S + J KfdV  . (21) 

Mere  the  factor  (1-v^/c)  can  be  placed  on  either  side  of  the  retarded  time 
derivative.  We  can  make  this  equation  look  more  like  a standard  Boltzmann 
equation  in  momentum  space  by  introducing  a modified  distribution  function 


F(T,p)  = (1  - -^)f  , f 


(1  - -f) 


Then  Equation  21  becomes 


where 


||  ♦ Vp  . (pF)  - S(T,p)  + /K*Fd3p’ 


1 

v dt 


fl  ‘if* 


K*  * K(p.p')  \ 


Now  from  Equation  1 wc  find  that  if  we  move  along  with  a particle. 


so  that 


dT  - (1  - -f)dt  , 


* 

*•3? 
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The  interpretation  of  K*  is  equally  simple.  The  scattering  kernel  K contains 
a collision  rate  and  an  angular  distribution.  Equation  25  shows  that  the 
collision  rate  in  retarded  time  is  increased  by  the  retarded-time  factor 
l/(l-v7c).  Particles  that  move  forward  with  close  to  c are  actc-d  upon 
by  forces  and  scattering  more  quickly  in  retarded  time,  because  the  real 
time  interval  is  longer  for  these  particles  than  the  retarded-time  interval. 
Note,  however  that  the  source  S does  not  acquire  the  retarded-time  factor. 

We  can  go  immediately  from  Equation  23  to  a completely  equivalent 
set  of  point  particles.  We  create  particles  according  to  the  probability 
distribution  in  S(T,p).  The  momenta  of  these  particles  change  at  the  rate 

= * usual  forces  , (28) 

1 - -5- 

c 

and  they  scatter  at  a rate  increased  by  a factor  l/(l-vz/c)  over  the  usual 
rate,  but  with  the  usual  angular  distribution.  In  adding  up  the  current 
densities  for  these  particles,  we  multiply  the  contribution  of  each  particle 
by  a factor  l/(l-vz/c)  because  that  factor  occurs  in  the  relation  (22) 
between  the  true  f and  the  modified  F.  This  is  the  method  used  in  CHAP, 
except  that  we  also  treat  the  scattering  approximately,  as  we  discuss  later. 

The  physical  explanation  of  the  retarded-time  factor  appearing  in 
the  current  contribution  of  each  particle  was  given  in  Reference  3.  It  is 
that,  for  two  Compton  electrons  bom  at  places  with  z-coordinates  differing 
by  AZq  and  having  identical  (but  displaced)  trajectories,  the  actual 
distance  apart  of  these  two  particles  at  any  given  time  is  AZ  = AZ^l-v^/c) 
as  indicated  in  Figure  2.  The  density  of  electrons  is  therefore  greater 


than  the  density  of  births  by  the  retarded-time  factor. 


Figure  2.  Explanation  of  retarded-time  factor  appearing  in  the  current 
contribution  of  each  particle.  Relation  between  AZ  and  aZq 
holds  in  the  limit  as  AZq  -*•  0. 


We  hope  this  discussion  of  the  retarded  time  equations  will  clear 
up  any  doubts  for  both  those  who  prefer  the  Boltzmann  approach  and  others 
who  prefer  the  particle  approach. 


18 


SECTION  5 

INITIAL  VALUE  OF  COMPTON  CURRENT 


For  a delta-function  pulse  of  gamma  rays,  we  can  calculate  exactly 
the  initial  value  of  the  radial  Compton  current.  Also  for  the  transverse 
Compton  current,  which  starts  from  zero  at  T = 0,  we  can  calculate  exactly 
the  initial  rise  rate.  These  two  quantities  are  not  affected  by  either 
energy  loss  or  scattering. 


The  contribution  of  a given  electron  to  the  radial  Compton  current 
is  proportional  to 


5.1 


z 


= Y(1+Y)X  • 


(29) 


Let  us  imagine  a magnetic  field  B in  the  y-direction  (see  Figure 
1),  which  will  lead  to  a current  in  the  x-direction.  The  contribution  of 
a given  electron  to  is  proportional  to 


6J 


x 


vx/c 


(30) 


The  average  of  this  expression  over  the  Klein-Nishina  distribution  vanishes 
by  synaetry.  The  derivative  with  respect  to  retarded  time  is 


Px  K 

(e-P2)2  ' 


(31) 


The  result  of  the  resistive  (energy  loss)  force  and  of  scattering  would  be 
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to  keep  the  distribution  symmetrical  in  the  angle  -p  about  the  axis  formed 
by  the  original  gamma  ray  direction.  Therefore  they  Jo  not  contribute  to 

Q 

the  average  of  W . Only  the  magnetic  force  will  contribute.  (we  assume 
here  that  there  is  no  electric  field,  although  we  could  also  calculate  its 
effect.)  for  the  magnetic  force. 


eB 

m 1 - vz/c 


;B  Pz 


™ t-  - P, 


VxA 


eB 

m 1 - v /c 

2 


eB  Px 


m c - p 


(32J 


We  then  calculate 


6J 


eB 1 

■'U- 


:-Pz)2  (c-Pz)3 


(53) 


Now  when  averaging  over  the  angle  <p. 


av(p^)  = av(pj)  = i av(P2+Py) 


1 2 2 
= 2 av(£  -l-pz)  • 

Thus  when  averaged  over  we  obtain 


(34) 


eB  1 1 (C  Pz5 
OJ  = T"  7 

X Cl  l , .3 

(e-pz) 


■ f T [(‘OX)3  - <>KX>] 


(35) 


Q 

To  find  the  initial  value  of  Jg  and  J , we  have  to  average 
Equations  29  and  35  over  the  Klein-Nishina  distribution  in  Xt  Equations  10 
and  11.  Tins  averaging  is  straight  forward,  if  a little  tedious.  The 
results  are,  in  MKS  units: 
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(36) 


•Jz(0)  = - NQecT2(Y)/T1  (f)  , 

jx(Oi  - ■»o“(?)T,tl,)/T1(Y)  • 07) 

Here  is  the  total  number  of  Compton  electrons  produced  per  unit  volume, 
Tj  is  given  by  liquation  13,  and 

T = (H-Y)h  * ~~ y - ?,n(U2Y)l  , (38) 

L (1+2y)~  Y“  Y J 


(1»2Y) ^ ,1  ♦ Y , 1 1 ♦ 2y  ? 

3 Y(1+2y)  2 „ 3 e'n<1'  *4 


For  comparison  with  later  numerical  calculations  we  record  here  values  for 

gamma  rays  of  energy  1.6  MeV  and  a magnetic  field  of  0.6  Causs  = 6 * 10-S 
2 

Weber/m  : 

- -Jz(0)/(N0ec)  = 7.70  , (40) 

Jx(0)/(N0ec)  = 26.54  = 2.80  x 108/sec  . (41) 


To  obtain  the  results  (36)  and  (37)  in  cgs  Gaussian  units,  replace 
ec  by  e,  and  eB/m  by  eB/mc. 


For  a given  flux  of  gammas,  Nq  is  itself  proportional  to  Tj , so 
that  Jz  and  are  proportional  to  T and  T . Graphs  of  these  quantities 
versus  gamma  energy  are  given  in  Figure  3. 
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SECTION  6 
ENERGY  LOSS 


l-'or  energy  loss  by  the  Compton  electrons  we  use  Bethe's  formula 
(Reference  9).  The  mean  change  of  energy  dW  per  track  length  ds  is 


■ ~2  d?  = 2™Zr0  T to- 
me p 

(42) 

Here  N is  the  density  of  atoms  of  atomic  number  2,  rQ  is  again 
electron  radius,  and 

the  classical 

[.]  - >.«  - s . yu2 

♦VioV- 

e 

(43) 

In  this  equation,  I is  the  mean  excitation  potential,  given  in 

for  various  elements.  For  air  and  aluminum: 

Reference  9 

air:  Z = 7.2  , 1 = 80.5  eV  , 

At:  Z = 13  , I = 150  eV  . 

(44) 

To  save  computational  time,  we  approximate  [l]  by 

[1]  2tn(~)  ♦ 3.42tn(p)  *1.71  , 

(45) 

« 15,80  ♦ 3.42tn(p)  for  air  , 

(46) 

*«  14.57  ♦ 3.42tn(p)  for  At  . 

(47) 

23 


This  approximation  is  accurate  to  i percent  tor  electron  kinetic  energies 
between  2s  keV  and  5 MeV. 


In  the  CHAP  method  we  imagine  that  a steady  force  equal  to  dh/ds 
acts  on  the  Compton  electron,  in  the  direction  opposite  to  its  velocity. 
Thu:,  th*  energy  of  the  electron  decreases  gradually  at  the  correct  average- 
rate.  ihis  treatment  ignores  the  fact  that  the  energy  loss  occurs  in  steps 
of  fluctuating  magnitude.  The  probability  distribution  of  energy  losses  w 
is  given  approximately  by 

2w 

Pfwjdw  ---■  -y--— — , (48 ) 

w0  ♦ w" 

where  w is  of  the  order  of  10  eV.  The  average  energy  Joss  per  inelastic 
collision  is 


w 


2 0 . , * . 

“V 

(i()  eV  for  K = 1 MeV 


(40) 


Thus  the  average  step  is  very  small  compared  with  the  energies  K of  the 
Compton  electrons.  However,  since  the  integral  of  wP(w)  gives  a logarithm, 
roughly  equal  amounts  of  energy  arc  lost  in  each  decade  in  w;  e.g. , about 
1/5  of  the  energy  of  a 1 MeV  electron  is  lost  in  collisions  that  lose  energy 
between  50  and  500  keV.  Thus  large  energy  losses  are  not  totally  negligible, 
fce  shall  see,  however,  that  at  the  center  of  the  high  altitude  RMP  source- 


region  (altitude  30  km),  energy  loss  has  only  a small  effect  on  the  peak 
Compton  current.  It  is  unlikely  that  fluctuations  in  energy  loss  could 
cause  changes  in  the  peak  current  of  more  than  a few  percent . 


In  CHAP,  ion  pairs  are  created  at  the  rate  of  one  ion  pair  per  85 
eV  lost  by  the  Compton  electron.  Secondary  ionization  then  proceeds  at  the 
rate  given  in  Reference  5,  until  there  is  one  ion  pair  per  34  eV  lost  by 
the  Compton  electron. 


SECTION  7 

COULOMB  SCATTERING 


The  differential  cross  section  of  a nucleus  of  .charge  Ze  for 
scattering  an  electron  into  angular  interval  dO  at  angle  0 from  its  original 
direction  is  (Reference  !)) 


, ,„2  2 r~  JnsiiiOdO 

d„  s r()  - . 

p q 

Here  q is  proportional  to  the  momentum  change, 


ISO) 


M = 2s  i n I ) (I)  u q --  2) 

- i.i  for  i)  « 1 . (Si) 

Note  that  the  differential  solid  angle  can  be  expressed  in  terms  of  q, 

dq~  - 2qdq  - 4s i n (~-)cos (y)dO  = 2sin0d0  . (S2) 


Thus 


do 


- 2.2 
4„z2r2  5bL 

' L *0  4 4 

P o 


(53) 


The  total  cross  section,  obtained  by  integrating  over  q,  is  infinite,  since 
even  distant  collisions  lead  to  some  scattering  for  a pure  Coulomb  potential. 
However,  screening  of  the  nucleus  by  the  atomic  electrons  reduces  the  scat- 
tering for  large  impact  parameters.  A good  way  to  take  into  account  the 
effect  of  the  screening  is  to  make  the  replacement 
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1 


\ 2 • (54) 

q a + q 

where  a is  effectively  the  minimum  angle  of  scattering.  Moliere  (Reference 
10)  found  a good  fit  to  the  scattering  from  the  Thomas-Fermi  atom  by  adding 
three  terms  of  this  form  with  different  a's  and  different  coefficients  re- 
placing unity  in  the  numerator.  We  have  picked  a single  value  of  a which 
gives  the  same  result  as  Moliere's  formulae  for  the  following  problem  (which 
will  be  used  in  the  next  section). 

We  wish  to  find  the  integral  over  the  differential  cross  section  of 
the  quantity 


1 - cosG  = 2sin2(|)  = ~ q2  , 


This  integral  is 


/<- 


cos0}do  = At\Z  r 


2^if 

0 4 2 /. 

p J I 


P4  2J  CaV)2 
0 


4"z‘ro  T I [*•<»  * ■?>  - — S’ 

p L • >‘VJ 


In  the  cases  of  interest  to  us,  a /4  will  be  very  small  compared  with  unity, 
so  that  we  can  approximate 

i[*n(1  + -t)  - — S]  - *«<§)  - \ = Hr1  • t57) 

a i ♦ V 

4 

Using  Moliere's  formulae,  one  can  again  do  the  integral  of  (l-cos0),  with 
considerably  more  work.  The  result  is 


j [l-cos8)do  « 4nZ2r2  ^ • 
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2 

(The  number  102  inside  the  logarithm  here  replaces  137  = ftc/e  in  less 
accurate  calculations  of  Coulomb  scattering.)  Comparing  the  Mol.iere  result 
with  liquations  56  and  57,  we  see  that  they  will  agree  if  we  choose 

,1/3 

= ITF  • <59> 

We  shall  use  the  Moliere  result  (58)  directly  in  the  obliquity 
factor  method  drived  in  the  next  section.  We  have  also  used  the  differential 
cross  section  (53),  with  the  replacement  (54)  and  with  a given  by  (59),  to 
construct  a Monte  Carlo  code  for  the  purpose  of  testing  the  accuracy  of  the 
obliquity  factor  method.  That  code  will  be  described  in  another  report. 

2 

One  often  sees  formulae  like  liquations  53  and  58  with  Z replaced 

by  Z(Z+1),  for  the  alleged  purpose  of  including  the  effect  of  scattering  by 

the  atomic  electrons.  This  procedure  may  be  approximately  correct  for  the 

larger  angle  scattering  (although  it  neglects  the  reduced-mass  correction), 

but  it  can  hardly  be  correct  for  the  small  angle  scattering,  where  the 

atomic  electrons  are  already  taken  into  account  in  the  screening.  Possibly 

some  form  like  Z(Z  + might  be  appropriate.  Since  the  correct  procedure 

* 2 

is  apparently  unknown,  we  leave  the  factor  V in  place. 


SECTION  8 

THE  OBLIQUITY  FACTOR 


In  Section  4 wc  derived  the  modified  Boltzmann  Equation  23  in 
retarded  time  and  momentum  space.  This  equation  conserves  particles  in 
momentum  space  (whereas  Equation  21  does  not)  and  is  directly  equivalent  to 
a set  of  particles  whose  accelerations  and  collision  rates  are  modified  by 
the  retarded  time  factor.  This  equivalence  is  exact  if  the  scattering 
operator  is  regarded  as  stochastic.  We  now  derive  an  approximate,  non- 
stochastic way  of  handling  the  scattering  of  the  Compton  electrons. 

We  consider  first  the  case  in  which  there  are  no  forces  (and  no 
energy  loss)  and  only  scattering  is  acting.  Then  Equation  23  becomes 

I F(P)  ' A(p.p')  dV  . (60)  ; 

J z'  I 

T.f  we  start  with  a particle  with  momentum  pQ,  then  F is  initially  a delta 
function 

F0(P)  ■ 'S(p-Pq)  . 

As  time  goes  on,  F(p)  will  spread  out  in  angle  around  p0,  but  with  no  change 

in  the  magnitude  p * pQ.  The  initial  angular  spread  will  be  small  because 

the  scattering  is  predominantly  small-angle.  Therefore  the  factor  (l-v’/c)  i 

-»  z 

will  vary  only  little  over  the  distribution  F(p’)  at  early  times,  and  may  j 

be  replaced  by  its  average  value  over  that  distribution.  We  thus  write  I 


JL 

3T 


Ffp) 


1 - Vz/c  /' 


K(p,p' )F(p' )d  p' 


(61) 


where  the  average  is  a function  of  T but  is  independent  of  p and  p' . 


In  this  approximation  the  central  angle  of  the  distribution  does 
not  change  from  the  angle  of  p(}.  Let  us  calculate  the  rate  of  change  of  the 
quantity  p defined  by 

p = J"cosaF(p)d3p  , (62) 


where  a is  the  angle  between  p and  pQ.  Multiplying  liquation  61  by  coso,  and 
integrating  over  p,  we  find 


dp 

dT 


1 - 


v2/c 


;// 


cosaK(p,p' ) F(p' )d3p'd3p 


(63) 


In  the  integral  over  p,  which  we  do  first,  we  can  choose  p'  as  the  axis 
of  spherical  coordinates.  Then 


cosa  - cosGcosS  - sinGsin8cos4>  , 


(64) 


where  the  angles  are  defined  by  Figure  4.  Since  K is  a function  only  of  the 
scattering  angle  6 and  not  of  the  second  term  on  the  right  in  Equation 
64  gives  no  contribution,  and  we  have 


(65) 


Figure  4.  Angles  for  Integration  of  Equation  S3. 
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In  these  variables,  the  second  integral  is  independent  of  p'  and  the  first 
integral  is  p.  Thus 


dp 

dr 


1 - 


v /c 
z 


/ 


cosGK(9)2TTsin0d0 


(66) 


Now  the  scattering  operator  K removes  particles  from  6=0  (i.e.,  it  contains 
a term  -6(6)),  and  puts  them  at  other  angles  6 > 0.  Since  it  conserves 
particles,  we  have 


/k(9)  2irsin0d0  = 0 . 


We  can  therefore  write  Equation  66  as 


dp 

dT 


J" (\-cosQ)  K(0)2itsin6d6  . 


(67) 


(68) 


Since  l-cos0  vanishes  at  0 = 0,  the  delta-function  part  of  K(Q)  gives  no 
contribution  to  this  integral,  and  K can  be  replaced  by  the  differential 
scattering  cross  section  of  Section  7,  multiplied  by  Nv  to  give  a scattering 
rate.  Wc  thus  obtain,  from  Equation  58,  and  noting  that 


v = c 


(69) 


the  result, 


1 _d_ 

c dT 


Unp 


- rrW  <”nz2'o  \ • (70) 


A group  of  particles  starting  out  with  velocity 


p 

v0  * c e * 


(71) 


will  have  an  average  velocity,  as  a result  of  scattering, 

V « v0p  . (72) 

This  reduction  in  velocity  affects  the  current  densities  and  the  relation 
between  real  time  t and  retarded  time  T for  this  group  of  particles.  It 
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does  not  affect  the  ionization  rate  directly,  which  depends  only  on  the 
magnitude  of  the  momentum. 


Let  us  now  turn  on  the  electromagnetic  and  resistive  forces.  Now 
the  magnetic  force  rotates  all  of  the  momenta  in  the  slightly -spread  distribu- 
tion (resulting  from  scattering)  about  the  magnetic  field  axis;  this  changes 
the  central  angle  of  the  distribution,  but  does  not  affect  the  angular 
spread  p.  The  resistive  force  reduces  the  magnitudes  of  all  the  momenta, 
but  also  does  not  affect  p directly.  The  electric  force  changes  the 
components  of  all  the  momenta  in  the  direction  of  the  electric  field;  this 
changes  both  the  central  momentum  and  the  angular  spread,  as  explained  in 
References  3 and  4. 


We  thus  arrive  at  the  obliquity  factor  method.  It  assigns  a 
central  momentum  p to  each  particle,  which  is  the  momentum  it  would  have 
in  the  absence  of  scattering.  Associated  with  p in  the  usual  way  are  the 
total  energy  e and  the  central  velocity  v.  The  equation  of  motion  in 
retarded  time  is 


it 

dT 


x usual  forces  . 


(73) 


1 'IT  ** 


The  equation  for  p is  Equation  70  (when  E = 0)  with 


vz  = vzp  . 


(74) 


The  contribution  of  a particle  to  the  current  density  is  proportional  to 


I* 


JLH. 


v e - p u 

1-fu 


(75) 


Since  p starts  from  unity  at  T » 0*  it  can  be  seen  that  scattering  does  not 

affect  the  initial  value  of  J , and  that  it  does  not  affect  the  initial  value 

z 

of  dJ^/dT  (for  a magnetic  field  in  the  y -direct  ion). 
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To  the  extent  that  scattering  is  dominated  by  small  angle  scatter- 
ing, it  is  clear  that  the  obliquity  factor  method  should  be  asymptotically 
correct  at  early  times,  i.e.,  it  should  give  the  first  order  effect  of 
scattering  correctly.  High  accuracy  at  late  times  is  not  expected  a priori , 
but  we  shall  see  what  we  get  in  the  next  section. 

In  earlier  discussions  of  the  obliquity  factor  method  we  used, 
instead  of  p, 

H = l/u  , £nn  = - fcnji  . (76) 

The  equation  for  n is  therefore 

“ - n * RHS  , (77) 

where  RHS  is  the  negative  of  the  right-hand  side  of  Equation  70.  In  our 
early  work  we  dropped  the  factor  n on  the  right  in  Equation  77,  on  the 
grounds  that  the  model  is  valid  only  when  n is  not  far  from  unity.  This 
method,  without  the  factor  n,  is  called  the  "old  method,"  whereas  Equations 
70  and  77  are  called  the  "new  method."  We  shall  compare  their  relative 
accuracies  in  the  next  section. 


SECTION  9 

COMPARISON  OF  RESULTS 


To  test  the  accuracy  of  the  CHAP  method  of  computing  Compton 
current,  we  shall  compare  its  results  with  those  from  Monte  Carlo  calcula- 
tions. In  another  report  we  shall  describe  two  types  of  Monte  Carlo  calcula- 
tions of  different  complexity.  The  simpler  calculation  treats  all  scattering 
as  small-angle,  using  a Gaussian  angular  distribution  of  scattered  particles; 
the  width  of  the  Gaussian  is  determined  from  the  total  scattering.  The  more 
sophisticated  calculation  divides  the  scattering  into  a small -angle  part 
and  another  part  not  limited  to  small  angles.  The  small  angle  part  is  as- 
signed a Gaussian  distribution,  whereas  the  large  angle  part  is  given  the 
screened  Coulomb  distribution.  The  sophisticated  calculation  is  considerably 
more  time  consuming.  We  have  used  both  methods  to  calculate  the  transmission 
of  monoenerget ic  electrons  through  aluminum  foils  of  various  thicknesses, 
and  have  compared  the  results  with  the  experimental  data  of  Marshall  and 
Ward  (Reference  11).  The  results  are  shown  in  Figure  5.  The  two  Monte 
Carlo  methods  give  ranges,  for  a given  transmitted  fraction,  which  differ 
by  not  more  than  about  5 percent.  The  experimental  results  are  very  close 
to  the  Monte  Carlo  results  for  small  foil  thickness,  but  show  ranges  up  to 
10  percent  larger  for  large  thickness  (low  transmissions).  We  do  not  know 
what  the  absolute  accuracy  of  the  experimental  data  is,  as  Reference  11  gives 
no  assessment  of  probable  error.  Spread  in  energy  of  the  incident  electrons 
would  make  the  tails  of  the  experimental  curves  extend  to  longer  ranges.  In 
addition,  crystaline  effects  in  the  aluminum  foils  are  not  accounted  for  in 
the  theory  of  multiple  scattering,  which  assumes  that  the  scattering  atoms 
are  randomly  placed.  This  effect  would  not  be  present  in  air.  Altogether, 
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curve  nearest  solid  curve  is  from  the  more 
cases.  Squares,  calculated  extreme  ranges 


the  agreement  is  quite  good,  and  we  shall  assume  that  the  sophisticated 
Monte  Carlo  represents  the  correct  effects  of  multiple  scattering  and  energy 
loss.  Since  the  simple  Monte  Carlo  gives  results  differing  by  only  a few 
percent  and  is  considerably  faster,  we  have  used  it  to  compare  with  the  CHAP 
method. 


All  of  the  calculations  were  done  with  gamma  rays  of  energy  1.6 
MeV  and  for  a transverse  magnetic  field  of  0.6  Gauss.  No  EMP  fields  (self 
consistency)  were  included.  Two  altitudes  were  used: 

-5  3 

altitude  = 30  km,  air  density  = 1.84  x io  gm/cm  ; 
altitude  = 20  km,  air  density  = 8.89  x io-5  gm/cnf*. 

Figure  6 shows  transverse  currents  at  30  km  altitude.  The  curve 
labeled  VAC  was  computed  by  the  CHAP  equations  but  with  the  resistive  force 
and  scattering  set  equal  to  zero,  i.e.,  it  represents  the  case  of  zero  air 
density.  For  the  curve  EL,  the  resistive  force  was  turned  on,  but  scattering 
was  omitted.  The  curve  labeled  CHAP  includes  both  effects.  We  see  that 
scattering  causes  a larger  effect  than  energy  loss.  The  curve  labeled  MC 
is  the  Monte  Carlo  result.  It  is  a few  percent  larger  than  the  CHAP  result 
at  times  of  several  nanoseconds.  We  shall  see  later  (Figure  9)  that  the 
ionization  rate  is  also  a little  larger  from  the  Monte  Carlo  calculations; 
these  two  errors  tend  to  cancel  in  determining  the  peak  electric  field,  which 
is  proportional  to  Jx/a  (o  is  the  conductivity).  We  see  that  all  the  curves 
approach  the  theoretical  slope  at  early  times,  although  scattering  causes 
departure  quite  early  in  retarded  time.  Hie  time  step  used  in  the  calcula- 
tions was  0.3  nanoseconds.  For  the  most  energetic  electrons,  the  real  time 
step  is  about  27  times  longer.  We  see  that  the  difference  between  the  old 
and  new  obliquity  factor  treatments  is  small  over  the  time  frame  presented, 
with  the  old  method  being  a little  closer  to  the  Monte  Carlo  results. 
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density;  EL,  energy  loss  Included,  but  no  scattering;  MC,  Monte  Carlo; 
CHAP,  new  obliquity  factor  method.  Circles  are  points  from  old  method 


Figure  7 shows  the  transverse  current  at  20  km  altitude.  Here  the 
Monte  Carlo  and  ('HAP  results  are  very  close  together  at  times  of  interest. 

_ g 

The  old  method  is  substantially  too  large  after  2 ' 10  second,  as  expected 
since  it  underestimates  the  scattering  when  ii  is  small  for  n is  large).  How- 
ever, this  error  would  not  affect  the  peak  EMP. 

The  transverse  currents  at  both  altitudes  are  shown  in  a linear  plot 
in  Figure  8. 

Figure  9 gives  the  ionization  rates  at  the  two  altitudes,  based  on 
instantaneous  production  of  one  ion  pair  per  34  eV  lost  by  the  Compton 
electron.  The  CHAP  code  takes  account  of  the  time  lag  for  secondary  ioniza- 
tion, but  we  have  omitted  this  lag  here  to  better  compare  the  calculations. 
(Including  the  lag  would  bring  the  Monte  Carlo  and  CHAP  results  a little 
closer  together,  since  both  curves  start  from  the  same  initial  value.)  At 
30  km  altitude,  the  Monte  Carlo  result  is  larger  than  the  CHAP  result  by  a 
little  more  than  the  error  in  the  transverse  current.  The  computed  peak 
EMP  from  CHAP  will  therefore  be  a few  percent  too  large. 

Figure  10  shows  the  radial  current  at  both  altitudes.  The  difference 
between  the  Monte  Carlo  and  CHAP  results  is  maximum  for  this  component. 
Fortunately,  the  radial  current  has  practically  no  effect  on  the  FMP  except 
at  points  very  close  to  the  burst,  where  EMP  is  usually  not  of  primary  con* 
cem.  The  radial  currents  all  start  from  the  theoretical  initial  value- 

As  a final  check  on  the  accuracy  of  our  calculations,  we  compare 
in  Figure  il  our  CHAP  results  with  some  older  Monte  Carlo  calculations  of 
Knutson  and  Morgan  (Reference  12).  These  calculations  were  made  for  20-km 
altitude,  but  with  slightly  different  and  B from  those  we  have  been 
using.  CHAP  was  run  with  these  revised  values.  Figure  II  shows  quite  good 
agreement  between  CHAP  and  the  Knutson- Morgan  Monte  Carlo.  Comparing  Figure 
11  with  Figure  8 (for  20-km  altitude),  we  see  that  the  relation  of  the  Monte 
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are  points  from  old  method. 
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Figure  9.  Ionization  rate  for  EY  = 1.6  MeV,  Bv  = 0.6  Gauss.  Curve  labels 
Indicate:  INITIAL,  initial  value;  JMC,  Monte  Carlo;  CHAP, 
results  from  new  obliquity  factor  method.  Circles  are  points 
from  old  method. 


Figure  10.  Radial  current  for  Ev  = 1.6  MeV,  By  = 0.6  Gauss.  Curve  labels 
indicate:  THEORY,  theoretical  initial  value;  MC,  Monte  Carlo; 
CHAP,  new  obliquity  factor  method.  Circles  are  points  from  old 
method. 


Carlo  curves  to  the  CHAP  curves  are  nearly  identical  in  the  two  figures. 
This  indicates  that  our  Monte  Carlo  results  and  those  of  Knutson  and  Morgan 
are  very  nearly  identical,  and  supports  the  reliability  of  all  of  the  cal- 
c"l at  ions . 


SECTION  10 

COMPARISON  WITH  SOLLFREY'S  RESULTS 


in  Figure  12  we  compare  Soil  Trey's  computed  transverse  currents 
with  our  Monte  Carlo  results  at  altitudes  of  20  and  30  km.  We  see  that 
there  are  substantial  differences.  Comparing  with  Figures  6 and  7,  we  see 
that  the  CHAP  results  are  much  closer  to  the  Monte  Carlo  results  than  are 
Sollfrey's  results.  This  is  disappointing  since  we  had  hoped  that  Sollfrey's 
calculations  would  provide  an  accurate,  independent  check  on  the  Compton 
currents. 


Not  having  gone  through  Sollfrey's  calculations  in  detail  (they 
are  quite  lengthy),  we  can  neither  confirm  them  nor  point  to  any  errors. 

We  do  raise  the  question,  however,  as  to  whether  the  series,  in  terms  of 
which  his  result  is  expressed,  is  convergent  or  only  semiconvergent . In 
this  connection,  we  point  out  that  his  first  term  is  closer  to  our  Monte 

Carlo  results  than  is  his  sum  JQ  + ♦ J2>  and  note  that  in  his  Figures  9a, 

b,  c,  the  series  does  not  appear  to  be  converging,  for  J2  is  generally 
larger  than  Jj.  It  appears  that  if  one  more  term  were  added,  the  result 
would  be  larger  than  the  vacuum  current  and  the  theoretical  initial  slope 
at  early  times.  In  our  view,  these  points  raise  serious  doubt  as  to  the 
accuracy  of  Sollfrey's  results.  We  would  like  to  se > these  questions  re- 
solved. 


Figure  12.  Transverse  current  for  Ev  * 1.6  MeV,  By  = 0.6  Gauss.  Curve  labels 
Indicate:  THEORY,  theoretical  initial  slope;  VAC,  zero  air  density; 
MC,  Monte  Carlo;  S,  Sollfrey's  computed  results. 


SECTION  11 
CONCLUSION 


We  have  developed  the  theory  of  the  Compton  current  in  the  presence 
of  the  geomagnetic  field.  We  have  derived  analytically  the  initial  value  of 
the  radial  Compton  current  and  the  initial  rate  of  rise  of  the  transverse 
Compton  current  (which  starts  from  zero  initially).  Neither  of  these 
values  is  affected  by  energy  loss  or  scattering  of  the  Compton  electrons, 
and  they  serve  as  checks  on  numerical  calculations.  We  have  explained  the 
approximate  but  fast  method  used  in  the  CHAP  code  for  calculating  the 
Compton  current,  including  the  effects  of  energy  loss  and  scattering.  We 
have  devised  an  accurate  Monte  Carlo  calculation  for  the  Compton  current  and 
have  shown  that  it  gives  good  agreement  with  experimental  data  on  the  trans- 
mission of  electrons  through  aluminum  foils.  We  have  compared  results  from 
CHAP  with  those  from  our  Monte  Carlo,  and  have  shown  that  CHAP  results  are 
within  a few  percent  from  the  Monte  Carlo  results.  Since  CHAP  errors  in 
Compton  current  and  ionization  rate  are  in  the  same  direction  and  about  the 
same  magnitude,  the  peak  Tsaturated)  electric  field  calculated  by  CHAP  should 
be  within  2 or  3 percent  of  the  correct  values.  We  have  shown  that  our  Monte 
Carlo  results  are  nearly  identical  with  similar  results  obtained  by  Knutson 
and  Morgan  in  one  case  available  to  us. 

We  have  compared  Sollfrey's  numerically  computed  Compton  currents 
with  our  Monte  Carlo  results,  and  found  discrepancies  as  large  as  50  percent 
at  important  times.  We  have  suggested  that  Sollfrey's  series  may  not  be 
convergent. 

Sections  1 through  9 of  this  report  mo  suitable  for  inclusion  in  a 
textbook  on  high-altitude  BMP. 
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